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Abstract 

We summarize a recent work done on the title's subject. First, we present the asymp- 
totic scheme of post- Newtonian (PN) approximation for general relativity in the harmonic 
gauge. Then, we discuss the definition of the mass centers and the derivation of equations 
for their motion, following that scheme. Finally, we briefly analyze the reason why a new 
term has thus been found in the equations of motion. 

1 Introduction 

The aim of this contribution ^ is to summarize a recent work iT on the approximate equations 
of motion relevant to general-relativistic celestial mechanics. The latter equations belong to 
the heart of the observational test of general relativity (GR), insofar as the first thing which one 
expects from a theory of gravity is to provide an accurate celestial mechanics. The first part 
of the paper will be concerned with the construction of the post-Newtonian approximation 
(PNA). Contrasting the "standard" scheme of PNA with the "asymptotic" one, which is 
used in this work, we will then present the application of the latter scheme to GR in the 
harmonic gauge and show some of the corresponding PN expansions. In the second part of 
the paper, the definition of the mass centers will be discussed and the way to derive explicit 
IPN equations of motion for the mass centers will be summarized. This derivation uses a 
recently-developed asymptotic framework for well-separated celestial bodies j2]. In addition, 
it is needed to exploit the fact that the main bodies of the solar system are nearly spherical 
and have approximately a rigid motion. We end the paper by showing the main new term 
found and discussing the reason why it occurs in this approach and did not in previous works. 



2 Post-Newtonian approximation 

The purpose of the PNA is to get approximate solutions of the Einstein equations for a weakly- 
gravitating system. In a further step, some integration should provide tractable equations of 
motion of the mass centers of the different bodies. Since a high accuracy is aimed at, the 
approximation scheme is an important point. The unknowns are: the (symmetric) metric 
tensor {gjj,u) (ten scalar unknowns), and the independent matter fields: the pressure p plus 
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the coordinate 3- velocity u = with x = (x*) the spatial position in the coordinate system 
(thus four scalar unknowns). We consider a barotropic perfect fluid, for which the proper 
rest-mass density p* and the internal energy density 11 are functions of the pressure p. The 
equations are the ten Einstein equations, plus the four equations of the gauge condition — 
here we adopt the harmonic gauge condition, which is usually used in relativistic celestial 
mechanics, and which is now indeed officially recommended [3] in this domain. 

2.1 Standard PN approximation 

This is the procedure followed by Fock in his classical book It is essentially equivalent 
to the one discussed later by Chandrasekhar [5]. It is still adopted in most of the current 
textbooks. In this approach, the metric {g^u), as well as the connection and the curvature 
tensor, are formally expanded in power series of c~^, where c is the velocity of light. However, 
the matter fields p,n, etc., are not expanded. As already noted by Futamase & Schutz jHj 
and by Rendall j2|, the latter point leads to a certain difficulty in interpreting the equations 
obtained. If one searches for asymptotic expansions in the usual sense, then all fields have to 
be expanded, moreover one has to properly define a small parameter e and to define a family 
(Se) of systems, in order that e can (conceptually) be made arbitrarily small. 

2.2 Asymptotic scheme of PN approximation 

This was first proposed by Futamase & Schutz |P, and this was for GR in the harmonic gauge. 
They derived the local equations (which are essentially the expansion of Einstein equations), 
though in a form which is not explicit enough for further use. As a matter of fact, "global" 
equations, i.e. ones for mass centers of extended bodies, were not derived by them, nor 
seemingly by anyone within this scheme, before the work In the present study, we start 
from Futamase-Schutz' family of initial conditions (slightly modified): 

pW(x)=eV'^(x), (1) 

p*(^)(x)=6V«(x), (2) 
u(^)(x) = eu«(x), (3) 



0. (4) 



[Here, g = det{gi^iu), and {g'^'^) is the inverse matrix of (g^Mu)-] Let be the system defined 
by this initial condition: we have indeed a family (S^) of systems. 

2.2.1 Expansion of the matter fields 

Let [M] and [T] be the mass and the time units for system Si. We change the mass and 
the time units for system S^, by adopting a large time unit [TJ^ = [T]/e and a very small 
mass unit [M],. = e^[M], where [M] and [T] are the units for system Si _8_. After this change, 
the initial data becomes independent of e, and we have just e = (we take c = 1 

in units for system Si). The initial- value problem still depends on e or more exactly on e^ , 
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since the Einstein equations involve only the square = e ^. This suggests stating Taylor 
expansions with respect to = c~^, starting from the zero-order term: 

p'=p'o+ + 0(c-4), u' = + u;c-2 + 0(c-^), (5) 

etc. In fixed units (the units of system Si), this is rewritten as 

P=\po+ Pie^ + 0(e')]e^ u = [uq + uie^ + 0(e^)]e, (6) 

etc. 

2.2.2 Expansion of the metric 

In a weak field, we expect that the metric g ^ rj, the flat metric. In particular, g should have 
an expansion with a first term of order zero (in fixed units): 

g= og + 0{e^). (7) 

Starting from the "relaxed" Einstein equations (valid in the harmonic gauge) allow us to 
deduce [T! that indeed 

og= r7 = diag(-l, 1,1,1). (8) 
With the time unit [T], = [T]/e, this is 

off' = diag(-e-2^ 1, 1, 1) = diag(-c2, 1, 1, 1). (9) 

To IPN order, one postulates Taylor expansions in c~^, up to and including the term, 
e.g. 

9m = -<? + 1500 + 2500 + (10) 

In the fixed starting units, this leads ^ to the expansion of the metric written by Weinberg 
[HI , including the g^i 's in . 

2.2.3 Expanded field equations: general form 

The expansion coefficients: PQ,Pi, i^oo; 2500) etc., are functions of the position x and of 
the "dynamical time" t' = et. They are by definition independent of e. Starting from the 
expansions of the matter fields and the metric, one deduces expansions for the other quantities, 
e.g. for the energy-momentum tensor T: 

T^- = (oT^- + iT'^- e2 ^ ^(^4)^ ^ ^ 2, n,o = 3, n,, = 4. (11) 

Again, the coefficients, e.g. qT^^ and iT^^ in (jlH) . do not depend on the parameter e. 
Therefore, if one enters these expansions into the Einstein equations and into the dynamical 
equations {T[l^ = 0), one gets each of these equations split to two separate equations, accord- 
ing to the order in e^: this is merely coefficient identification in a polynomial function. 

The general form of the split equations, i.e., with the coefficients qT^'^, iT^^ , not with 
their explicit expressions, is the same as for Weinberg. This is because he too expands T. 
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2.2.4 Explicit expansions for a perfect fluid 

We enter the expansions of the metric and the matter fields, into the expression of T for a 
perfect fluid: 

T^'^ = {^* + p)U^'U'' + pg^"" (c = 1) (12) 

{U^ = dx^/ds is the 4-velocity, ^* = p*{l + 11) is the proper energy density). This leads to 
the explicit expansion of tensor T. E.g.: 

oroo = /.5, iroo = /ii + /.5(ug-2ci.). (13) 

($ is the Newtonian potential.) Weinberg ,9^ does not expand the matter fields /i*,u, etc.. 
He denotes by p what we note p* in Eq. ()12() . but then he writes, for instance, 

or°° = p, ir°° = p(u2 - 2$). (14) 

One may tentatively interpret the unexpanded matter fields of Eq. (|14j) within the asymptotic 
scheme, as being the first- approximation fields. This interpretation means, for instance: 

/5SW = ^(1) = (^0 + /^!^^)^^' usw = u(i) = e(uo + uie^). (15) 

Then the coefficients o^sW' i^sw CSJ) do depend on the small parameter e. (This is 

called a "composite expansion.") But then, it is difficult to see why the field equations could 
be split. On the other hand, if one does not split the equations, while expanding some of the 
independent variables (here the metric is expanded), then one has not enough equations. More 
precisely, at the IPN order, the expansion of the metric component goo uses two independent 
unknowns, Eq. ()1U() . Hence, for this component, a problem arises (IJ if the corresponding 
Einstein equation is not split to two equations. 

2.2.5 Explicit expanded equations for a perfect fluid 

Inserting the asymptotic expansion of T, Eq. (|13|) and the like, into the dynamical equation 
T-fif^ = 0, gives the equations of fluid dynamics. The zero-order equations turn out to be those 
of classical fluid dynamics: continuity equation and Euler's equation. They are exact for 
the zero-order fields! That continuity equation expresses the mass conservation at the order 
zero, not at PN order. However, the order-one time equation, combined with the zero-order 
equations, allows one to check that mass is conserved also at the IPN approximation. (In 
fact, mass is exactly conserved for an isentropic perfect fluid in GR.) Due to the expansion of 
matter fields, the order-one equations are definitely different from the corresponding equations 
of Weinberg [0] , which coincide with those of Chandrasekhar j^j . 

3 The definition of the mass centers and their equations of 
motion 

3.1 Definition of the mass centers: motivation 

In any relativistic gravity, any form of energy must 

• contribute to the gravitational field; 

• be subjected to its action. 
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Therefore, a question arises: which energy density should one choose as "the" correct weight 
function (f) so as to define relevant mass centers? 

The energy density (j) adopted as the weight function should 

• i) be well-correlated with the luminous density; 

• ii) ensure that the velocity of a body's mass center is the average velocity of this body's 
constituents. 

Condition ii) is satisfied if and only if (j) satisfies the continuity equation: 

dtcf) + div(0u) = 0. (16) 

Therefore, the rest-mass density p (in the global frame), which satisfies these two conditions 
and which seems to be the only energy density to do so, should be chosen as the weight 
function (j). A more standard choice [e.g. Will ^0], Eqs. (6.21) and (6.25)] is to consider a 
density p' = p plus: i) the density of internal energy, ii) the energy density associated with 
the self Newtonian potential of the relevant body, and iii) the density of the kinetic energy 
associated with its motion with respect to a local frame attached to the mass center of that 
body. (However, several well-known works also take "(/> = p", e.g. Fock Brumberg jl5j.) 
Since p' — p is of order p x GMa/ {(?ra) where is the size of body (a), the two different 
choices can lead only to small and, most importantly, nonsecular differences. 



3.2 Formal definition of the mass centers 

The masses and the mass centers are defined with the global rest-mass density p: 



Ma= / pdV, Maa= / pxdV, (17) 

where Da(t) is the spatial domain occupied by body (a) (a = 1,...,A^) in the considered 
(harmonic) coordinate system. 

At the IPN approximation, p is approximated by p(i) = po + Pic~^, hence 



M(i) = mO + ^ [ po dV, Ml ^ j pi dV, 



(18) 



with 



Mii)a(i) = / p(i)X dV = M% + M>ic-2, (19) 

JDa 

M°ao= [ po^dV, Miai= [ p^x dV. (20) 

JDa JDa 

Note that and ao are the Newtonian mass and mass center, because the zero-order fields 
obey the Newtonian equations. The definitions (|19|) and (|20j) seem natural. Also note, how- 
ever, that the IPN position of the mass center, a(;^), is not equal to ao + aic^^. In other 
words, ai is not the IPN correction to the position of the mass center. 
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3.3 IPN equations of motion of the mass centers: general form 



To get the PN equations of motion of the mass centers, the local PN equations of motion 
(obtained as described in S I2.2.5|) are integrated in the domain D^. More precisely, integration 
of the space components of the local PN equations of motion gives the sought equations for 
motion of the PN mass center a^^) of body (a). At the order zero, we get the Newtonian 
equation of motion: 

M°ao = - / pV<^>^''Uv (21) 

(<!>(") is the external part, for body (a), of the Newtonian potential For the IPN correction, 
we get: 

_ -P + J" + - Afjao 
^(1) - ao - , (22) 

where I", J",K" are integrals of PN fields over D^, hence are structure- dependent. Thus, a 
priori, the equations of motion are so. 



3.4 Good separation between celestial bodies 

For many self-gravitating systems, in particular for our solar system, the main bodies are well 
separated. This means that one may introduce a small parameter thus: 

7]Q = max — <C 1 (rb = radius of body domain D;,). (23) 

aj^b |a — b| 

To account for this, we assume that the IPN equations are accurate enough for the system 
of interest, S. In other words, S is replaced by the corresponding IPN system S'. (This is in 
order to avoid worrying with two small parameters, e and ry.) Then, we build a family of IPN 
systems, (S'^), with S"'" = S'. This is defined by a family of initial conditions which ensure 
that 

• (r^i,)^ = |aQ — bpi = ord(r/^-^) for a b, 

• (ao)'' = ord(r/i/2)^ 

• (n("))^ = ord(r/V2)^ n('i) = rotation velocity . 

(The zero-order velocity field is assumed to be a rigid rotation, see just below.) 



3.5 Justification of assuming rigidly-rotating bodies 

Even at Newtonian order, rigid rotation cannot be exact, due to the time-varying influence 
of the other bodies, which produces tides. However, if one writes Euler's equation for a 
rotating body in a well-separated system, one finds that, neglecting 0{r]^), the body's internal 
equilibrium is not affected by the other bodies. This accuracy for the internal equilibrium 
turns out to be enough to get equations of motion up to rj^ included. To this accuracy, it is 
hence consistent to assume rigidly-rotating bodies. 
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3.6 Quasi-spherical bodies 

In addition to having nearly rigid rotation, the big celestial bodies are nearly spherical, but 
not exactly so: cf. the case of Jupiter. This can be described in an asymptotic framework by 
coupling the good separation with the quasi-sphericity: we assume that 

= 0{rf) (a = l,...,iV; i,A; = 1,2,3), (24) 

where the 7}^^^ 's are the eigenvalues of the inertia tensor of body (a). 
As a result, it is shown that the rotation velocity rate verifies 

d^"^ = Oirj% (25) 
which may be neglected for the ord(r/^) equations of motion. 

3.7 Equations of motion (well-separated, rigidly-rotating, quasi-spherical 
bodies) 

In this framework, integrals I", J", K*^, entering the equation for IPN correction to the accel- 
eration of the mass center, Eq. ((HJ, can be computed up to rj^ included. Inserting the result 
in this equation gives the explicit equations of motion in the asymptotic scheme. The New- 
tonian and IPN equations are thus separated. To compare with the standard Lorentz-Droste 
|llj (or Einstein-Infeld-Hoffmann ,12]) equations, one needs to group both together. This 
turns out to be feasible. As a result, we find ^ that the acceleration of the mass center of the 
generic body (o) is exactly the Lorentz-Droste acceleration A]^^, plus a IPN correction 
that cancels if all bodies are exactly spherically symmetric (this term is likely to be extremely 
small in the solar system), plus a IPN correction depending on the rotation velocity of the 
body and on its internal structure: 

.ao + A- + 0(r?4) + 0(c-4). (26) 

Here, ja is the spherical inertia moment, which does depend on the density profile, hence 
on the internal structure; and {lu^"'^)'^ = Q^^g^Q^j^j^ /2 is the square of the angular velocity. 

[fi^^^ri^") is the spatial tensor (in Cartesian coordinates) with ik component Q^j^Q^^^ . Note 
that ao is the Newtonian acceleration ()21() of the body considered.] We note that this new 
term is order 3 in tj. It is clearly a self-acceleration term. This new term does not seem to be 
negligible for the giant planets, and, the author conjectures, it will lead to secular effects. 

3.8 Why is the new term here and was absent in previous works? 

There are essentially two different approaches to the equations of motion. 

1. There are works based on the local integration of the IPN field equations obtained 
in the global reference frame for a perfect-fluid system, as is based the present work. The 
previous works that belonged to this category {e.g. Fock j^, Misner et al. Spyrou [Tl] . 
Will [im, Brumberg [IS]): 



(a) (a) 



a(i) 



LD 



7a 
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• used the standard (Fock-Chandrasekhar) approximation scheme summarized in Sub- 
section 12.11 The present author beUeves to have shown that this scheme is not really 
compatible with an asymptotic interpretation ^ Ej. Of course, both the standard 
and the asymptotic PNA are just approximations, and thus, in a sense, they are both 
"wrong". However, it seems to this author that the asymptotic scheme has serious 
chances to approximate accurately the exact behaviour — because, precisely, it is fully 
consistent with asymptotic analysis; 

• did not involve either an asymptotic description of the good separation between bodies. 
(This point is discussed in Ref. 0.) The same lack led, in the first version 17^ of the 
equations of motion of the scalar theory investigated by the author, to neglect terms 
which were later found to be numerically significant. 



2. More recently, there have been works based on multipole series expansions with local 
frames (Damour, Soffel, and Xu jl8| I19j. Racine and Flanagan [20]). These works have led 
to tractable equations of motion only for two models: the "monopole model", which gives 
the Lorentz-Droste acceleration, and the "monopole-dipole model" , which adds a (negligible) 
spin-dependent term. Here we note that these authors use a rather complex kind of multipole 
series expansion, so that it is difficult to get an asymptotic estimate of the omitted remainder. 
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